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The general theory of disclinations developed earlier is applied to the special case of a straight 
disclination line. First the geometrical fields are found, such as the defect loop densities which cor- 
respond to Mura's new concepts of "plastic distortion" and "plastic rotation," the basic plastic fields 
(strain and bend-twist), the defect densities (dislocation and disclination), the characteristic vectors 
(Burgers and Frank), and the incompatibility. Then the static fields are found for the isotropic case, 
such as the displacement, total distortion, basic elastic fields, and the stress. It is shown that the dis- 
clination axis is moved by adding a dislocation to the disclination line. All these special results for the 
straight disclination line are shown to satisfy the general equations of the theory. As corollaries the 
following topics are also treated: 

1. The finite and infinitesimal straight disclination dipole. which can be biaxial or uniaxial. It 
resembles the straight dislocation line. 

2. The dislocation models of the straight disclination line and of the finite disclination dipole. They 
are terminating dislocation walls (tilt and twist). 

3. The compensated disclination line and the bent dislocation wall. 

4. Finally we show analytically a special case of a dislocation ending on a disclination. 

Key words: Burgers vector; continuum mechanics; defect; dipole; disclination; dislocation; distortion; 
Green's tensor; incompatibility; loop; plasticity; strain; Volterra. 
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1 . Introduction 

In previous papers [1,2, 3] 1 we have developed the general linear theory of stationary discli- 
nations in an anisotropic or isotropic body. In the present paper we wish to specialize some of these 
results to the case of discrete straight disclination lines. 

In 1907 Volterra [4J gave the elastic fields of straight hollow core dislocations and disclinations 
running along the axis of a long isotropic cylinder. He called these defects "distortions." Later Love 
[5] ventured to call them "dislocations." This name has been widely accepted, but its meaning has 
been narrowed to translational dislocations, because of the early importance of this type of defect 
in plastic deformation. Therefore Frank [6] introduced the new name "disinclination," later changed 
to disclination, to describe the rotational dislocation. Furthermore the word "distortion" has since 
come to stand for the gradient of the displacement, and it would therefore lead to confusion if we 
still tried to ascribe to it Volterra's original meaning as the combination of dislocations and dis- 
clinations. Still it would be useful to have a name for this combination, and we have adopted 
the word defect for that purpose. In this terminology, then, the present paper deals with discrete 
straight defect lines. This word usage does not exclude the consideration of "point defects," 
because they can be regarded as combinations of discrete dislocation loops [3]. 

We would here also like to point out that it has become customary to refer to "disclination 
theory" when really "defect theory" is meant (i.e., the theory of dislocations and disclinations). 
We have done this in our paper titles to emphasize the fact that the new results in defect theory are 
obtained because disclinations have been introduced. 

Furthermore there are at least three independent ways of defining the two types of defects, 
i.e., in terms of (1) continuous distributions, (2) discrete lines, or (3) infinitesimal loops, which lead 
to essentially different descriptions of the disclination. This means, for example, that a discrete 
disclination line contains a certain amount of dislocation density, or a certain different amount of 
dislocation loop density. In this paper we shall identify these quantities for the case of straight 
disclinations. 

In section 2 we recall that Volterra also summarized the characteristic terms in the displace- 
ment fields for straight defects. These results have to be slightly modified to obtain equilibrium 
everywhere. 

Section 3 reviews straight dislocations to show the approach used in specializing the general 
theory. It also serves as a reference for the later sections. 

In section 4 we apply the general results of [2] to our special case of discrete straight disclina- 
tions. We find the defect loop densities which are Mura's "plastic distortion" and "plastic rotation." 
Then we identify the basic plastic fields from which all other results follow. We find the defect 
densities and show that they satisfy the necessary general relations. The Burgers and Frank vectors 
are examined, and we conclude with relations for the incompatibility tensor. 

Section 5 treats the isotropic statics of straight disclinations. From the given displacement 
we derive the basic elastic fields, and show that they satisfy the field equations with the defect 
densities. Then we find the stress and show that it maintains equilibrium. 



1 Figures in brackets indicate the literature references at the end of this paper. 
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In section 6 we find how the fields are modified when the disclination axis is moved away from 
the origin. We show that it is equivalent to adding a dislocation to the disclination. 

Section 7 treats the infinitesimal straight disclination dipole line. 

In section 8 we give the dislocation model of the straight disclination. We show that it is a 
terminating dislocation wall (in general tilt and twist). It has the same displacement, elastic strain 
and stress as the corresponding disclination line. In terms of the dislocation model we find more 
precise interpretations of Mura's new concepts. 

Section 9 treats the finite straight disclination dipole line. We show that it reduces to the 
infinitesimal case of section 7 when the width approaches zero. 

In section 10 we derive the results for the straight compensated disclination line, using two 
different defect surfaces. 

In section 11 we apply the results of section 10 to the bent dislocation wall, which is made up 
of two semi-infinite dislocation walls meeting at right angles. 

Section 12 shows analytically how a dislocation can end on a disclination. 

Throughout the paper we assume that the body is infinitely extended. The geometric results 
are valid for the anisotropic case, but we have obtained the static fields only for the isotropic body. 

We shall freely refer to the results of the three previous papers [1, 2, 3], which will be denoted 
in the present paper as I, II, and III. 

2. Volterra's Equations 

Weingarten's theorem [1, 2] states: On following around an irreducible circuit in a multiply- 
connected body, the displacement and rotation change by an amount that represents a rigid body 
motion, if and only if the classical elastic compatibility conditions are satisfied throughout the body. 

This theorem led Volterra to conceive of the defects he called "distortions." They have sub- 
sequently also been called the most general Weingarten-Volterra dislocations. Volterra then 
proceeded to calculate the properties of such straight defects lying along the z axis. He first found 
the terms in the displacement fields that characterize such defects. Nabarro [7] considered these 
so important to dislocation theory that he quoted them directly in his book (eq (2.20)): 

2ttu x = (b x — dzxZ+d xy y) tarn' 1 (y/x) + i (- b y - d yz z — ad xy x)\n (x 2 + y 2 ) , 

27Tu y = (bij — d X yX+d yz z) tan _1 (y/*) + \ (b x — d zx z — ocd xy y) In (x 2 + y 1 ) , (2.1) 

2ttu z = (b z — d yz y+d zx x) tan -1 (y/x) + i (d yz x-\- d zx y) In {x 2 + y 2 ) , 

where a is an elastic constant. These displacements maintain equilibrium under no body forces 
except at the z axis. It will be recognized that the terms with coefficients b represent possible dis- 
placements for dislocations. The terms with coefficients d will be seen to be possible displacements 
for straight discrete disclinations. 

These equations have to be modified to obtain equilibrium under no body forces also at the 
z axis. This can be done in several ways. One method is to calculate the body force from (2.1) 
by the usual equations of elasticity. It is found to be singular at the z axis. Then the displacement 
field of an equal but opposite force along the z axis is derived, and added to (2.1) to get the desired 
result. We do not show the details in this paper. Another method is to use the general expressions 
from III. Even though straight defects do not satisfy the necessary boundary conditions used there, 
it is still possible to find the correct functional dependence of the fields from the general expres- 
sions. The details are given in appendix C. 

Once a set of displacement functions is found "by hook or crook," it is easy to verify 
that they are a solution by showing that they maintain equilibrium everywhere, including at singular 
points, and that they include the characteristic Volterra terms (2.1). Equilibrium is checked simply 
by deriving the stress and showing that it satisfies the equilibrium equation. 
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The development in this paper can also be put in terms of the following general question: 
Given the displacement as an empirical function of the coordinates representing experimental 
data, what is the system of defects that would produce such a displacement? To answer this 
question we have to determine how to divide the resulting fields into elastic and plastic parts. For 
a discrete defect this will be simple to do, because the most logical candidates for the basic plastic 
fields are the terms that are singular on a particular surface, as we shall see for the dislocation in 
section 3.2 and for the disclination in section 5. However, by another division we can also find a 
different system of defects that will give the same displacement, as we shall also see for the "dis- 
location model" of the straight disclination in section 8. 

3. Review of the Straight Dislocation 

In this section we apply the general results of section II 3.2* to the special case of a straight 
dislocation. 

3.1. Geometry 

The plastic distortion for a constant displacement across a surface S and the dislocation 
density for the corresponding discrete dislocation along the line L are given by the general expres- 
sions, (113.21, 113.23),* 

/8£(r)=-8*(S)6i, (3.1) 

a p i(r) =8 p (L)bi, (3.2) 

where L is the closed boundary of S, 8k(S) and 8 I} (L) are the Dirac delta functions on the surface 
S and the curve L, respectively, defined by (All-12) in appendix A, and bi is constant, identified 
as the Burgers vector of the discrete dislocation line. For the special case of this paper we wish 
to take I as a straight line and S as a semi-infinite plane bounded by L. In particular, let L run 
along the x 3 axis and S be the part of the X\X% plane for negative X\. To satisfy the right-hand rule 
we take the normal of S in the negative %2 direction. These relationships are illustrated in figure 1. 
Then the only nonvanishing components of (3.1 - 2) become 

/8£ = 6iff(-*i)8(*0, (3.3) 

a 3l = bi8(x l )8(x 2 ) 9 (3.4) 

according to (A15-16). The plastic distortion is singular on the plane S, and the dislocation den- 
sity is concentrated on the X3 axis. The three different cases for / = 1, 2, 3 are illustrated in figures 
2-4. In figures 2 and 4 the plane S is also the slip plane, because it contains the Burgers vector. 
For consistency we see that these results satisfy the definition, (15.2), (113. 4), 

ol iA = - e pmk ji p kl m , (3.5) 

by (A6). We see also that (3.4) satisfies the continuity condition, (15.4), (113.24), 

otpi tP = 0. (3.6) 

Now take for a Burgers circuit X a circle of radius a in the X\X>z plane centered at the origin (fig. 1). 



*Section II 3.2 refers to section 3.2 in reference [2], and the symbol (113.21) refers to eq (3.21) in reference [2]. 
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Figure 1. Geometry of the discrete straight defect line in 
the Cartesian coordinate system XiX-zX.i. 

The defect line L is along the x* axis and the defect surface S is the semi-infinite 
part of the JC1JC3 plane for which x\ is a negative, with a normal in the negative x 2 
direction. For the Burgers circuit A, which links L t and the Burgers surface o~, 
which is pierced by L, we have chosen a circle of radius a in the Xi*j plane centered 
at the origin. 




Figure 2. Edge dislocation described by the dislocation 
density a 3 i = 6i8(jci)8(x 2 ) . 





FIGURE 3. Edge dislocation described by the dislocation 
density 0:32 = WUi )8(x 2 ). 



FIGURE 4. Screw dislocation described by the dislocation 
density a 3i = b :i 8{x\ )8{x 2 ). 
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Then we have by (3.3) and (A23) 



&p p kl dL k = -hb,H(-x,)S(x 2 )dL 2 = b, 



(3.7) 



in agreement with (113.22), confirming that bi is the Burgers vector. We also have by (3.4) and (A22) 

I a p idS p = \ bid (xi) 5 (x 2 )dS3 = 6/, (3.8) 

J (T J (T 

in agreement with (113.25). So in this section we have found the special forms of the plastic distor- 
tion and dislocation density (3.3^1) for a straight dislocation running along the x 3 axis, and these 
equations are in complete agreement with all the general results for geometry. 

We conclude this section with expressions for the incompatibility tensor. It is in general given 
by, (15.3), 

t)pq = — (CqnKXpl, n)(pq). (3.9) 

For our special case the only nonvanishing components are 

Vi3=-ib 3 8( Xl )8'(x 2 ), (3.10a) 

T)23=ifc 3 o'(*i)o(* 2 ), (3.10b) 

T)33 = 6i8(*,)8' (x z ) -b 2 8' (xi)8(* 2 ) , (3.10c) 

from (3.4). We note that they satisfy the continuity condition, (14.2), 

Vpq,p=0. (3.11) 

The results in this section apply to an infinitely extended body, whether isotropic or anisotropic. 

3.2. Isotropic Statics 

In this section we find the static fields for a straight dislocation in an infinitely extended iso- 
tropic body. The displacement is given by (appendix C) 



f-»'[£ + i^] + C(!^)[< l -*'>"' + 2] 



2tt' 



(3.12a) 



(3.12b) 



(3.12c) 



where v is Poisson's ratio. These expressions are well-known for a dislocation line of Burger 
vector bi running along the x 3 axis. We have used a cartesian coordinate system ^i^ 2 ^3 , as shown in 
figure 1. The cylindrical coordinates p and </) are defined in appendix B and shown in figure 5. By 
restricting the angle </> to the range (—77, 7r) it becomes discontinuous on the semi-infinite plane 
S shown in figure 1, which then becomes the natural plastic displacement surface of section 3.1. 
Note that we have labeled (3.12) as the total displacement. In the literature these expressions are 
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Figure 5. Relation between Cartesian coordinates xi 
X2 and cylindrical coordinates p and (j). 



and 



The angle </> has the range (— tt, tt), assumes the values shown at the ordinates, 
and has a jump of 2tt across the negative x x axis. 



usually regarded as the elastic displacement which is then taken to be multiple-valued on account 
of the angle (/). In the framework of the present theory, however, it would not be correct to regard 
(3.12) as the elastic displacement. For, if the elastic displacement is prescribed, then there can 
be no dislocations (c.f. 12.1, 12.2, and 15.6, c.f. also footnote 10 in II). But if the total displacement is 
given the dislocations are still unspecified. We see from (3.5) that they are specified when we de- 
termine how the total distortion is split into an elastic and a plastic part. Intuitively we can see 
that it is unreasonable to split the displacement into an elastic and a plastic part in the present 
special case because of the fact that the angle </> cannot be split into a continuous and a discrete 
part. In general we shall contend that if a total field quantity of a straight defect contains the 
angle (/>, it is unreasonable to split it into an elastic and plastic part. This idea will also be applied 
in section 5. The results (3.12) show a logarithmic dependence for the displacement at large dis- 
tances from the dislocation line, in agreement with table 113. 

The total distortion is defined in terms of the total displacement and consists of an elastic and 
plastic part, (12.1), (II3.9), 



Pli = ul k = p k , + Pi 



(3.13) 



On differentiating the displacement (3.12), the terms containing </> will produce expressions singular 
on S, by (B4). These expressions correspond exactly to the plastic distortion (3.3) found in the pre- 
vious section. This apparent coincidence is of course achieved by our special choice of S and the 
range of 4>. These terms occur only because we chose the angle <£> to be discontinuous rather than 
multiple-valued. Alternatively, if we had no a priori knowledge of the plastic distortion, we could 
assign the singular terms to it. These terms intuitively represent a surface along which a plastic 
displacement has taken place, as discussed in section 3.1. So we find for the elastic distortion of 
the straight dislocation 
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/3 31 = 0, 



(3.14c) 
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4tt( 



l-v)[ 



(l-2^- + 2-- 2 
P 2 P 4 



62 



4tt(1 -v) . 



(3-2»)— -2 



c 1 ^ 2 " 1 



(3.14d) 



£22 = - 
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*•" 2wp«' 
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2irp 2 
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(3.14f) 

(3.14g) 
(3.14h) 

(3.141) 



using (B3-7). By contrast with the plastic distortion (3.3), we see that the elastic distortion consists 
of functions that are continuous everywhere except at the X3 axis. This is reasonable to expect for 
an elastic field. To check the consistency of our results, it is easily shown that the elastic distortion 
and the dislocations density (3.4) satisfy the field equation, (15.6), (113.10), 



€pmk fikl, m ~ OLpl , 



(3.15) 



by (B8-10). This result then confirms that (3.14) represents the elastic distortion of a discrete 
dislocation line running along the x-s axis. 

Next we want to show that equilibrium is maintained everywhere. We first find the elastic 
strain, which is in general given by, (113.13), 



eki — fi(ki). 



(3.16) 



Hence we have from (3.14) for the elastic strain of a straight dislocation 
6, 
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e 2 2 = — 
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The dilatation e = ekk is 
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At large distances from the dislocation line the elastic strain varies inversely with the distance, 
in agreement with table 113. Hooke's law for an isotropic body is in general given by, (III3.2), 



(T k i = 2G e k i + 



\-2v 
Hence we find for the stress of the straight dislocation 
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(3.19a) 
(3.19b) 
(3.19c) 
(3.19d) 
(3.19e) 
(3.19f) 



These expressions are also well-known in the literature. Finally, the equation of equilibrium are, 
(112. 1), 

fi = -cr k i, k , (3.20) 

where ft is the body force per unit volume. From (3.19) and (B8-10) we find 

fi = (3.21) 

everywhere, including at the x 3 axis where (3.19) is singular. So the above fields maintain equilib- 
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rium without body forces, including at the X3 axis. The displacement (3.12) contains the character- 
istic Volterra terms (2.1) for a dislocation. Therefore, as discussed in section 2, we may conclude 
that (3.12) correctly represents the displacement of a dislocation. 

We conclude this section with expressions for the elastic rotation of the dislocation, which is in 
general given by (III4.6), 

(3.22) 



(3.23a) 
(3.23b) 
(3.23c) 



In this section we have summarized the well-known expressions for the displacement and the 
elastic fields of a straight dislocation in an infinite isotropic body. 

4. Geometry of Straight Disclinations 

In this section we shall apply the general geometrical results of section II5.2 to the special 
case of a straight disclination. For a discrete defect along the line L, which has been formed by 
plastic displacement across the surface S, the dislocation and disclination loop densities are in 
general given by, (115.12-13), 

/3*(r)=-8,(S){6/+e/ 9r a,(x r -<)}, (4.1) 

<,(r)=-8*(S)ft 9 , (4.2) 

where bi is the dislocation Burgers vector, fl q the Frank vector, x? a point through which the 
disclination axis passes, and 8 k (S) is the Dirac delta function on the surface S, defined by (A12) in 
appendix A. The quantities (4.1-2) were introduced by Mura [9], who called them the "plastic 
distortion" and "plastic rotation." For our special case of a straight disclination running along the 
X3 axis, we take S to be the semi-infinite ^1^3 plane for negative X\ with normal in the negative X2 
direction, as illustrated in figure 1, and set bi=0. Furthermore, for simplicity at this point, we 
let the axis of rotation pass through the origin, so that x® = 0. The effect of moving the axis will be 
investigated in section 6. Then the only nonvanishing components of (4.1-2) become 

^ l = €i Qi n Q XrH(-x l )8(x2), (4.3) 

ip*=a < fl(-x 1 )8(x 2 ) 9 (4.4) 

according to (A16). For convenience, we write these results out in full 

Pf^CLiXsHi-xOSixi), (4.3a) 

j8*= (n 3 x 1 -a l x 9 )H(-x 1 )8(x 2 ), (4.3b) 

P* 3 =-n 2 xiH(- Xl )8(x 2 ), (4.3c) 
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rf^ftiftC-xOSta), (4.4a) 

<p* = £UH(- ^1)8(^2), (4.4b) 

^*-H,//(-x I )8(x 2 ), (4.4c) 

where we have used (A8). The basic plastic fields (plastic strain and bend-twist) are in general 
given by, (115.14-15), 

K^ q = Uklqfikl, m+<Pm- (4.6) 

For our special case we easily find the only nonvanishing components for the plastic strain 

ef 2 =i02*3#(-si)8(s 2 ), (4.7a) 

e% 2 = (Cl3Xi-£liX 3 )H(-xi)8(x 2 ), (4.7b) 

e£ 3 = — ifi2XiH(— Xi)8(x2), (4.7c) 

from (4.3), and by (4.6) we find that the components of the plastic bend-twist are 

icfj— m«ff(— *i)8(jfc) f (4.8a) 

K^=-Ul2X 1 H(-xi)8 , (x2)^il 1 H(-xi)8(x 2 ) 9 (4.8b) 

= 0, (4.8c) 

(4.8d) 

kZ 2 = Q,2H(-Xi)8(x 2 ), (4.8e) 

< 2 = 0. (4.8f) 

*f 3 = ^2^38(^1)8(^2), (4.8g) 

K£ 3 =-&l 2 *3#(-*i)8'(*2)+ft3# (~*i)8(* 2 ), (4.8h) 

4 i --in 2 //(-^ 1 )8Cx 2 ), (4.8i) 

from (4.3-4), where we have used (A6). The defect loop densities (4.3-4) and the basic plastic 
fields (4.7-8) are singular on the plane S, figure 1. This is the plane on which a discrete amount of 
plastic displacement has taken place to create the disclination line along the Xz axis. 

The general expressions for the dislocation and disclination densities are given by, (115.19-20), 

a P t(r) = 8 p (L){bi+ei qt a q (xr--x°)}, (4.9) 

pq (r) = 8 p (L)Cl q , (4.10) 

where 8 P (L) is the Dirac delta function on the curve L, defined by (All). For our special case, with 
bi = x$ = 0, the only nonvanishing components of the defect densities are 

OL-M— eiqr£lqX r §{X\)S(x2) , (4.11) 
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by (A15). Again, for convenience, we write these results out in full 

a M = 112X38(^)8(^2), 
a 32 = — flix 3 8 (xi ) 8 (xi ) , 



03i = a 1 8(x 1 )8(x 2 ) 9 
Ss2 = n 2 8(x l )d(x 2 ), 

e 3 3=a 3 s(x 1 )8( X2 ) 9 



(4.12) 

(4.11a) 
(4.11b) 
(4.11c) 
(4.12a) 
(4.12b) 
(4.12c) 



where we have used (A8). The defect densities are concentrated on the x 3 axis. We illustrate the 
three different types of disclinations for q=l, 2, 3 in figures 6-8. Figures 6-7 represent twist 
disclinations because the Frank vectors Hi and 2 are normal to the disclination line, whereas 
figure 8 represents a wedge disclination because the Frank vector fl 3 is parallel to the disclination 
line. From (4.12) we see that the twist disclinations correspond to the off-diagonal components of 
the disclination density, whereas the wedge disclination corresponds to the diagonal component 
(see table III). Note also that the discrete straight twist disclination contains a certain amount of 
dislocation density, given by (4.11a-b). The significance of these expressions can be interpreted 
as follows. Consider a slice of material parallel to the X\X 2 plane at the point #3. Locally near the 





FIGURE 6. Twist disclination with Frank vector Qi. 

It is described by the defect densities cx 32 — — fti* 3 8(.*i )8(x 2 ) and d 3 \ 
n 1 8(x l )8(x 2 ). 



FIGURE 7. Twist disclination with Frank vector fl 2 . 

It is described by the defect densities a3i = n 2 Jt 3 8(jci) 8 (x 2 ) and 6 32 = Ct 2 8(xi)d(x 2 ). 
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FIGURE 8. Wedge disclination with Frank vector ;{ . 

It is described by the disclination density 8 33 — il 3 8(x\ )8(x 2 ). 



X3 axis the displacement jump of the disclination cannot be distinguished from that of a dislocation, 
even though there is a rotation jump. In figure 6 the deformation and the displacement jump near 
the point x* is similar to that of the dislocation in figure 3 with the density given by (4.11b). Similarly 
in figure 7 the deformation near the disclination line is similar to that of the dislocation in figure 2 
with the density given by (4.11a). From this point of view it is therefore quite reasonable to regard 
the discrete disclination line as consisting of a mixture of dislocation and disclination density. 
This is a particular illustration of the general assertation we have already made in II. Note that 
the dislocation density implies a dislocation line with changing Burgers vector. 

By contrast Anthony [10] does not regard (4.11) as a dislocation density, but as an integral part 
of the disclination density. The advantage of this approach is that the disclination is unambiguously 
defined, whether it is described in terms of continuous distribution of defects, or a discrete line. 
However, this approach becomes cumbersome when we try to describe a terminating dislocation 
as in section 12. Anthony has also treated several special cases of straight disclinations. For the 
wedge disclination his equations (27~28) and (30) correspond to our equations (4.12c) and (4.11c). 
We shall compare his results for the twist disclination with ours in section 6. 

For consistency it is straightforward to verify that the basic plastic fields (4.7~8) and the defect 
densities (4.11 — 12) satisfy the general definitions, (16.3, 16.1), (II4.1— 2), 



&pl — €pmk \€frf , m • ^-klqKmq ) -, 



(4.13) 



Opq— € PmkK kqm , 



(4.14) 



by (A6) and (A8). Furthermore, it is easy to verify the general relations (II6.27~28) between the 
loop densities (4.3~4) and the defect densities (4.11-12) 



OLpl — — € P mk (0w. m + *klq<P * ) , 



(4.15) 
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9pq=-Cpmk<P$ q , m , (4.16) 

also by (A6) and (A8). We also see that the defect densities (4.11-12) satisfy the general continuity 
conditions, (16.8, 16.2), (115.27-28), 

OLpi,p + ei pq d pq =0, (4.17) 

e pq , P =0. (4.18) 

We next wish to calculate the total Burgers vector which is in general defined in terms of the basic 
plastic fields by, (II4.5), 

Bl= ~J > X Kl- € l<irK P kQ Xr)dL k , (4.19) 

where A is the Burgers circuit, figure 1. For our special case we find from (4.7-8) by a somewhat 
tedious calculation that 

Bi = 0, (4.20) 

where we have used (A8-9). The general relation for a discrete defect line is, (115.17), 

B l = bi-ei qr ^ ( pc% (4.21) 

and with our special assumption that bi = x^. = this reduces to (4.20). The total Burgers vector Bi 
measures the displacement jump at the origin, which vanishes in figures 6-8. The total Burgers 
vector can in general also be expressed in terms of the defect loop densities, (116.31), or the defect 
densities, (II4.7), 

Bl = - ^ (j8fi ~ €l Q r^ Q X r ) dL k (4.22) 

= I (oipi — ei qr 0p q Xr)dSp, (4.23) 

where o" is the Burgers surface, figure 1. From (4.3-4) and (4.11-12) it easily follows that these 
equations lead to the same result as (4.20). Now we take for the Burgers circuit A. a circle of radius 
a in the x t x 2 plane centered at the origin (fig. 1). Then we have by (4.8) and (A19-21) 

-&K p kq dL k = tt q , (4.24) 

in agreement with (115.18), confirming that Cl q is the Frank vector. We also have by (4.4) and (A23) 

-j> <p% q dL k = -j> il q H(-x 1 )d(x 2 )dL 2 = n q , (4.25) 

in agreement with (116.32), and by (4.12) and (A22) 

J 6p q dS p = J n q 8(x 1 )8(x 2 )dS 3 = £l Q , (4.26) 

J a J a 

in agreement with (115.30). The jump in displacement across the surface S is in general given by, 
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(115.1,115.5), 

[a/(r)] =Bi + €i qr £l<pcr 

= bl + €lqr£lq(x r — X*) . 

In our special case with bi = x° r =0, or by (4.20), this reduces to 

[ui(r)] = ei Qf {l Q Xr. 



(4.27) 
(4.28) 

(4.29) 



This jump agrees with those depicted in figures 6-8. So, in this section we have found the special 
forms of the defect loop densities (4.3-4), the basic plastic fields (4.7-8), and the defect densities 
(4.11-12), for a straight disclination running along the x$ axis, and these equations are in complete 
agreement with all the general results for geometry from II. 

We conclude this section with expressions for the incompatability tensor. It is in general 
given by, (16.6), (114.12), 

Vpq=— (€qnlOLpl,n+0pq)(pq)- (4.30) 

For our special case we find the only nonvanishing components from (4.11-12) 

T) 13 =-fti8(*i)8(*2), (4 31a) 

7}23 = -il 2 8(Xi)8(X2) (4.31b) 

r)^ = a l x&8 , (xi)8(x2)^^kx 3 8(xi)8 f (x 2 )-il38(xi)8(x2), (4.31c) 

by (A9). We note that these results satisfy the continuity condition, (14.2), (114.13), 

7) M , P =0. (4.32) 

The results in this section apply to an infinitely extended body, whether isotropic or anisotropic. 



5. Isotropic Statics of the Straight Disclination 

In this section we find the elastostatic fields for a straight disclination in an infinitely extended 
isotropic body. If the Yolterra expressions (2.1) are properly modified to maintain equilibrium 
everywhere, then the displacement is given by (appendix C) 
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where v is Poisson's ratio. We have used a Cartesian coordinate system xi 9 x 2 , £3, as shown in 
figure 1, and the cylindrical coordinates p and </>, as defined in appendix B and shown in figure 2. 
This is the displacement field of a straight disclination line running along the x 3 axis, with Frank 
vector ft 9 . By restricting the angle (/> to the range (— 7T, tt) we shall find that the semi-infinite 
plane S becomes the natural plastic displacement surface of section 4. The asymptotic dependence 
of the displacement (5.1) goes as pin p, in agreement with table 113. 

The total distortion is in general defined in terms of the displacement by, (12.1), (II3.9), 



PL*"*. 



(5.2) 



In performing these differentiations in the displacement (5.1) we shall regard three-dimensional 
space to consist of the product of the two-dimensional x x x 2 plane and the independent one- 
dimensional x 3 line. For x x and x 2 we use the results for generalized derivatives of appendix B. 
Since x 3 appears only linearly in (5.1) the classical and generalized derivatives with respect to x 3 
coincide. So we find for the total distortion, using (B3-7), 
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(5.3d) 
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(5.3h) 
(5.3i) 
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Note that the total distortion contains terms which are singular on S and which correspond exactly 
to the dislocation loop density /3 A * given by (4.3). This result is in agreement with the general 
equation (115.34). Therefore we may be tempted to regard these terms as the plastic distortion 
(3%! and consequently find the elastic distortion from (3.13). This is the point of view taken by 
Mura [9]. We shall investigate the consequences of this approach in section 8. Strictly speaking 
we are at liberty to split the total distortion in any way we wish into an elastic and a plastic part. 
However, the present theory then implies that there are no disclinations. We note, though, that 
the total distortion (5.3) still contains the angle </>, in contrast to the dislocation case. Since </> is 
discontinuous at the surface S, it is not possible to split it into a part that is continuous everywhere 
except at the*.} axis, and a part that is localized only on the surface S. Therefore the total distortion 
cannot be split into an elastic part that should be continuous everywhere except at thex ;i axis and a 
plastic part that is localized on the surface S of plastic displacement. We shall find, however, that 
the basic total fields no longer contain the angle (/> and so it will be easy to find an intuitively reason- 
able way to split them into an elastic and plastic part. 

We first consider the total strain which is the symmetric part of the total distortion and con- 
sists of an elastic and plastic part, (114.14), 



e li = P(ki)= ekl + e W 



(5.4) 



If we use (5.3) in this expression we shall find that the terms which are singular on S correspond 
exactly to the plastic strain e p kl given by (4.7). So we find for the elastic strain of the straight 
disclination 
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(5.5g) 



The elastic strain consists of functions that are continuous everywhere except at the x* axis, which 
is reasonable to expect for an elastic field. These expressions also show a logarithmic dependence 
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for the elastic strain at large distances from the disclination line, in agreement with table 113. 
To find the total bend-twist we first find the total rotation which is in general given by, (1115.11), 



<* T q = UkiqPl r 



or, from (5.3), for our special case by 
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coT=^-ln 2Xl H(- Xl )8(x 2 ), 



Q 2 ip 



"* 2n ' 



(5.7a) 
(5.7b) 



ft 1*2*3 fl 2*1*3 l n Ul v s , x ^3<P 

znp 2 2-rrp 2 I 2tt 



(5.7c) 



These expressions still contain the angle $, which prevents us from splitting the total rotation 
into an elastic and plastic part in a reasonable way. The total bend-twist is in general defined in 
terms of the total rotation and can be split into an elastic and plastic part, (114.15), 
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(5.8) 



In carrying out this differentiation on the functions in (5.7) we shall find that the resulting terms 
singular on S correspond exactly to the plastic bend-twist k£ given by (4.8). So we find for the 
elastic bend-twist, using also (A6), (B3^4), and (B8), 
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The elastic bend-twist, like the elastic strain, no longer contains the angle </> and is continuous 
everywhere except at the X3 axis. Alternatively, if we had no a priori knowledge of the basic plastic 
fields (4. 7~8), we could assign to them the terms in the basic total fields that are singular on 5, 
because these surface singularities intuitively show where the plastic displacement has taken place. 
In this manner we would then also have been able to obtain the basic elastic fields (5.5) and (5.9), 
starting only with the expressions for the displacement (5.1), as discussed at the end of section 2. 
To check the consistency of our results, it can be verified that the basic elastic fields (5.5) and 
(5.9) and the defect densities (4.11-12) satisfy the field equations, (16.3, 16.1), (114.17-18), 



€ P mk ( Ckl, m H" GklqKmq) — OL p i t 
Gpmk K kq, m pq » 



(5.10) 
(5.11) 



by (B8-12). These results then confirm that (5.5) and (5.9) represent the basic elastic fields of a 
discrete disclination line running along the X3 axis. Furthermore it can be verified that the elastic 
strain (5.5) and the incompatibility tensor (4.31) satisfy the field equation, (114.22), 



" € pmk €qnl€ kl , m n 7) pq , 



(5.12) 



by(B8-12). 

Next we want to show that equilibrium is maintained everywhere. We first find the stress from 
Hookes' law (3.18) and the elastic strain (5.5) 
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(5.12b) 
(5.12c) 
(5.12d) 
(5.12e) 
(5.12f) 



If we substitute these results into the equations of equilibrium (3.20) we find by (B5-10) that 

/*=0, (5.13) 

everywhere, so that the above fields maintain equilibrium without body forces, including at the 
X3 axis. The displacement (5.1) contains the characteristic Volterra terms (2.1) for a disclination. 
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Therefore, as discussed in section 2, we may conclude that (5.1) correctly represents the displace- 
ment of a straight disclination line. 

In this section we have presented the static fields, i.e., the displacement, total distortion, basic 
elastic fields (strain and bend-twist), and stress, of a straight disclination line in an infinite isotropic 
body. 

6. Movement of the Disclination Axis 

If we do not restrict the disclination axis from passing through the origin as in sections 4-5 we 
obtain more general results. These can be found by removing the restriction x° r =0, which was used 
in those sections to find the explicit formulas from the general expressions of II. In this section 
we shall discuss a few selected expressions which are of particular interest. They will show that the 
effect of moving the disclination axis is equivalent to adding a dislocation, as was already pointed 
out in II. 

6.1. Geometry 

The general expression for the dislocation density of a discrete defect line is given by (4.9). We 
shall now consider the special case of a discrete disclination line lying along the x 3 axis, but whose 
rotation axis passes through the general point x?. Then the only nonvanishing components of the 
dislocation density become by (A15), setting 6/=0, 

a:i!=€i qi Ci Q (x r — x° r )8(x l )8(x 2 ), (6.1) 

instead of (4.11). On the other hand, the disclination density of our discrete disclination line with 
its axis shifted from the origin remains the same as (4.12). When written out in full we have for (6.1) 

a 3 i = [ft 2 (%3-* 3 °) +ft*S]8(«i)8(*2) , (6.1a) 

a 32 = [-Ct 3 x° l -Ct l (x 3 -x 3 )]d(x l )d(x 2 ), (6.1b) 

a 33 = [-fti^+ft2<|8(*i)8(*2), (6.1c) 

where we had used (A8). These results of course reduce to (4.11) when jt® = 0. In the present case 
we find that the total Burgers vector no longer vanishes, 

Bi = -ei qi £l&» r , (6.2) 

from (4.21) with fc/ = 0. So the displacement jump across the surface S for our special case becomes 

[ui(r)] = ei qr Uq(xr — x° r ), (6.3) 

from (4.27-28), or, when written out in detail 

[u l (r)]=tl 2 (x 3 -x»)-n 3 (x 2 -x< 2 ), (6.3a) 

[u 2 (r)]=a 3 (x 1 -x\) -Oi(*a-j$), (6.3b) 

[u 3 (r)]=Ct l (x 2 -x 2 )-Ct 2 (x l -x 1 ). (6.3c) 

The total Burgers vector (6.2) is the displacement jump at the origin. Anthony [10] has treated the 
special case of the twist disclination, Hi , for the case when the axis passes through the point 
x$ = x 2 (Z), %3 = 0. To find the results for his case we set il 2 = Cl 3 = in the above equations. The 
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only nonvanishing components of (4.12), (6.1), and (6.3) then are 

0si = Oi8(xi)8(x s ) 9 

a^ = — Ct l x 2 (Z)d{x l )8(x 2 ), 
[u 2 ]=-rL lX3 (Po), 

[u 3 ]=a 1 (x 2 (Po)-x 2 (Z)) 9 



(6.4) 

(6.5a) 
(6.5b) 
(6.6a) 
(6.6b) 



where we have set x 2 = x 2 (Po) , #3 = #3(^0). These results correspond to Anthony's equations (31), 
(33), and (34). Hence our results are in complete agreement with Anthony. 

If we compare (6.1) with (4.11) we see that the axis has been moved from the origin to the 
point x® by adding the dislocation density 



aM = — ei Qr ft q x i) t 8(xi)d(x 2 ) 



(6.7) 



to the disclination line of section 4. This is in agreement with (115.24). From (3.4) this corresponds 
to adding a dislocation line with Burgers vector 



bl = — eiqrTLqX® 

to the disclination line, in agreement with (115.23). In detail (6.8) is 



(6.8) 

(6.8a) 
(6.8b) 
(6.8c) 



So the disclination axis is moved from the origin to x® by simply adding a dislocation with Burgers 
vector (6.8) to it. 



6.2. Statics 

The fields of a straight discrete disclination along the^ 3 axis whose rotation axis passes through 
the point x° r can be obtained by adding the fields of the dislocation line in section 3.2 with the 
Burgers vector given by (6.8) to the fields of the disclination line in section 5. For example, the dis- 
placement is obtained by substituting (6.8) into (3.12) and adding it to (5.1) 
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<= ft 4^^-^a^" (lnp - 1) J- n 4^^ + MT^ x2(lnp - l) J- (6 - 9c) 

This is the displacement of a straight disclination line running along the X3 axis, with a rotation 
axis that passes thru the point x%. Other field quantities can be calculated in the same way, or can 
be derived from (6.9) by the same method as used in section 5. 

7. The Infinitesimal Disclination Dipole 

In this section we wish to apply the general expressions of section 117 to the special case of 
a straight disclination dipole line. The biaxial dipole is obtained by moving the whole basic straight 
disclination (line and axis) by an infinitesimal distance f m . The uniaxial dipole is obtained by 
moving only the line and keeping the axis fixed. In both cases the negative of the basic disclination 
is left at the X3 axis. 

7.1. The Biaxial Dipole 

In general the defect densities of a biaxial dipole are given by, (II7.1— 2), 

<Xpl= — €m<Xpl,m, (7.1) 

Opq = C, m Opq, m , (7.2) 

where ct p i and 6 pq are the defect densities of the corresponding basic defect. For the special case 
of a discrete straight biaxial disclination dipole line the basic defect densities are given by (4.11- 
12). So for the conjugate dipole the only nonvanishing components of (7.1—2) are 

a§,=-€nrO«[*r6i8 i «(p) + fr8(p)], (7.3) 

05,= -«*fm8, m (p), (7.4) 

where 8(p) is the two-dimensional Dirac delta function defined by (B9). These results could of 
course also be found from (117. 3- 4) and (A15). In these expressions the tensor ftqf m represents the 
strength of the dipole. When written out in full (7.3-4) become 

af 1 = -n 2 [^fi8'(xi)8(jC2)+«368Ui)8'(x2) + 68(^)8(^2)], (7.3a) 

c^=n 1 [*aft8^i)8(*»)+*rfi8(*i)8'(*«)+68(xi)8(ft)], (7.3b) 

% = -Oi[fi8'(*i)8(x,) + {,8(*i)8' (*,).], (7.4a) 

% = -n,[fi8'(*i)8(*,) + ft8(*08'(xi)], (7.4b) 

fla=-n3[fi8'(*i)8(*,) + fi8(xi)8' (*,)], (7.4c) 

where we have used (A8— 9). 

In general the displacement of the biaxial dipole is given by, (117. 5), 

a?=-£m< m , (7.5) 

where u J is the displacement of the corresponding basic defect. For our special case the distortion 
of the basic straight disclination line is given by (5.3), and so the displacement of the conjugate 
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dipole is 



a 2 6*3 



M i 1 — 7^ ^T (1—2^)— — 2 — — — - — — -\ (6 — 2v)— — 2 

1 47r(l — v) L P 2 P 4 J 47r(l — i/) L P 2 p 4 



+ 



(Oifi + Orfi)*, 



47t(1 — v) 






47r(l — i^) 



(7.6a) 



*?= 



i#=- 



£l\£\x* 



47r(l — ^) 

47r(l — i/) 



(3-2v)- 2 -2-^ 

P 1 P 4 



47r(l — v) [ p 2 p 4 



(1-2^) 



P 2 P 4 J 



^i** -| (l 2 ^-ft 3 f 2 



l-2i/) lnp+4 
P 2 



iiifi r *i 

"3=4^)l (1 - 2 " )ln ^ 



^ + _ a 2 6 



4w(l-v) 



4tt(1-^) 
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(7.6b) 



(7.6c) 



where we have ignored the terms in (5.3) which are concentrated on the surface S. We presume 
that the displacement is not a state quantity because it contains the angle </>. The displacement has 
a logarithmic dependence at large distances from the dipole line, in agreement with table 113. 

7.2. The Uniaxial Dipole 

To derive the relations for the uniaxial disclination dipole we first find the results of moving 
the axis of the basic disclination line by an infinitesimal distance fm, keeping the disclination line 
fixed. The dislocation density of this defect is given by, (II7.7), 



<*&= £ mda P J dx^, 



(7.7) 



whereas the disclination density vanishes, 6jj Q =0. For our special case of a straight disclination 
line we find from (6.1) that the only nonvanishing components of (7.7) are 



in agreement with (117. 9) and (A15). The displacement of the above defect is, (117.11), 

ui = t; m dulldx m , 
which, for our special case, becomes 



i*?=-(n 2 & 



-o*> [I 



+ • 



XiX 2 



2tt- 4tt(1-i')p 2 J 4tt( 
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10b) 
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(7.10c) 



from (6.9). If we compare (7.8) and (7.10) with (3.4) and (3. 12) we see that the above results are the 
dislocation density and displacement of a discrete straight dislocation line with a Burgers vector 
given by 



bl — — e\qi^lq^r, 



(7.11) 



in agreement with (117.13). 

We now find the relations for the discrete straight uniaxial disclination dipole line. In general 
they are given by the difference between those in section 7.1 and the above. For example, the dis- 
location density is found from, (117.14), 

*£ = *#-«£• (7.12) 

Hence from (7.3) and (7.8) we find that the only nonvanishing components of (7.12) are 

OL^ — — €lqrfl<pCri;m8,m(p), (7.13) 

in agreement with (117.15) and (A15). When written out in full this expression becomes 



afi=-n^,Ki8'(*i)8(x,) + f,8(*i)8 f (* f )]-nafrf(xi)8(x,), 
a& = n 1 * 8 [68'(*i)8(x«)+f,8(xi)8'(x,)]+nrfi8(* 1 )8(*,), 
a&=(ni£,-n,fr)8(*i)8(*«), 



(7.13a) 
(7.13b) 
(7.13c) 



where we have used (A8-9). The disclination density of the uniaxial dipole is the same as for the 
biaxial dipole, L = QP q , given by (7.4). The displacement of the uniaxial dipole is given by, (117.16), 



v n *~n- 



By (7.6) and (7.10) this becomes for our special case 
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p 2 J 
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+ (0,fe-Orfi) 
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(7.15c) 



We note that these expressions no longer contain the angle <p, and so u\\ is a state quantity. The 
asymptotic dependence at large distances is the same as the biaxial dipole, and shown in table 113. 
Furthermore the displacement of the uniaxial dipole u]\ does not contain any contribution 
from the wedge disclination, i.e., ft 3 is missing from (7.15). We may therefore conclude that it is 
possible to find the displacements of an edge dislocation from those of a biaxial wedge disclination 
dipole. Specifically if we set 



we find for (7.11) 



n, = ft 2 = 0, 



=n 3 £ 2 



-n 3 £i 



(7.16) 



(7.17a) 



(7.17b) 



If these results are substituted in (7.6 a-b) we find (3.12 a-b). This method of obtaining the edge 
dislocation displacements was first pointed out by Eshelby [11]. 

8. The "Dislocation Model" of the Straight Disclination 

8.1. Geometry 

As discussed in section 119 the dislocation model of a discrete defect line along the curve 
L is obtained by eliminating the disclination loop density (119. 1) 



<Pkq z 



0, 



(8.1) 



from the corresponding defect surface S that spans L. We shall illustrate this for our special case 
and so find the dislocation model of a straight disclination. There are of course many possible 
choices for the surface S, but we consider only the one used in the previous sections and illustrated 
in figure 1. 

The dislocation loop density remains the same as before, (4.3), and becomes the plastic dis- 
tortion, (II9.2). So it has the only nonvanishing components 



13!;, = /?* = € IgrilgXrH (-Xi)8( X2 ), 

or, when written out in full 

0&=j8* = (Osxi— Oix»)ff(-xi)8(* t ), 
j8£=/3*,=-ft2X,i/(-x 1 >8(x8). 

The dislocation density of the dislocation model is in general given by, (II9.7), 



(8.2) 

(8.2a) 
(8.2b) 
(8.2c) 
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a p i(r)=8p(L){bi+ei qr i),q(xr — x°)} + 8i(S)£l p — 8 p i8 k (S)Cl k - 



(8.3) 



For our special case, with bi = x^. = 9 this reduces to 

an = fl2#(-*i)8(* 2 ), 

a 3l = rL 2 X3&(xi)d(x 2 ) , 
ai2 = -Ct 1 H(—x 1 )8(x 2 ), 

«22 = 0, 



a 3 2 = — £lix 3 S (x 1 )8(x 2 )— n 3 //(— *i)8 U2), 



«13 = 0, 



«23 = 0, 



a 3 3 = n 2 ^(— *i)8(s 2 ), 



(8.4a) 
(8.4b) 
(8.4c) 
(8.4d) 
(8.4e) 
(8.4f) 
(8.4g) 
(8.4h) 
(8.4i) 



where we have used (A15-16) and (A8). The dislocation model consists of two parts, the same dis- 
location density along L as for the corresponding discrete disclination, given by (4.11), and 
a constant distribution over the surface 5, given by the terms with the Heaviside function. 
We have sketched the three different type of dislocation models for fti, H 2 , and fl 3 in figures 9-11. 
Rather than the constant dislocation distribution over S implied by (8.4), which is somewhat dif- 
ficult to draw, we have shown the discrete distribution that most closely approximates it, and is 
somewhat easier to visualize. These figures represent the dislocation models of the straight dis- 
clinations shown in figures 6-8. Figures 9 and 11 represent terminating tilt walls and figure 10 a 
terminating twist wall. The disclination density of course vanishes, (II9.8), 



e Da =o. 



(8.5) 




Figure 9. 



.••-.••'•••• 



Dislocation model of the twist disclination in 
figure 6. 




The tilt wall on the surface S (fig. 1), which terminates at the x 3 axis, is repre- 
sented by a set of parallel edge dislocations. The arrows are Burgers vectors and 
the attached short lines show the orientation of the extra half planes. This defect 
is described by the dislocation densities oc i2 = — fti//(— X\ )&{x 2 ) and oc 3 2 = — 
n t x»8(xi)8(xt). 



FIGURE 10. Dislocation model of the twist disclination in 
figure 7. 

The twist wall on the surface S, which terminates at the x 3 axis, is represented 
by a crossed grid of screw dislocations. The arrows are Burgers vectors and the 
attached short lines show the orientation of the extra half planes for the edge dis- 
locations. This defect is described by the dislocation densities a u ~a 33 = Q,2H 
(—Xi)8(x 2 ) and a 31 = il 2 x 3 8(xi )8(x 2 ). 
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FIGURE 11. Dislocation model of the wedge disclination 
in figure 8. 

The tilt wall on the surface S, which terminates at the x% axis, is represented by 
a set of parallel edge dislocations. The arrows are Burgers vectors and the attached 
short lines show the orientation of the extra half planes. This defect is described 
by the dislocation density a 32 = — 3 //(— Xi)8(jt 2 ). 



For consistency it is easily verified that (8.2) and (8.4) satisfy the definition (3.5). It is also 
easily checked that the continuity condition (3.6) is satisfied by (8.4), in agreement with (119. 9). 

We conclude this section with expressions for the incompatibility tensor. If we substitute 
(8.4) into (3.9) we find the only nonvanishing components 



7]i 3 = — Cii8(xi)d(x 2 ), 
T723 = -n 2 S(*i)8(% 2 ), 



(8.6a) 



(8.6b) 



7)33 = 0,1X38' (xi)8(x 2 ) + £l2X3d(xi)8' (x 2 ) —£l3&(xi)d(x 2 ). 



(8.6c) 



These results are identical with the incompatibility tensor of the corresponding disclination given 
by (4.31). Since the incompatibility is the source of elastic strain, (114.22) and (114.31), this means 
that the dislocation model has the same elastic strain as the corresponding discrete disclination 
line. 

8.2. Statics 

According to (119.11-12) the total displacement and distortion of the dislocation model are 
identical to those for the corresponding disclination line. They are therefore given by (5.1) and (5.3). 

To determine the elastic fields it is necessary to know how to split up the total fields. The dif- 
ference between a discrete disclination line and its dislocation model simply results from splitting 
the total fields up in a different way. For a discrete disclination we have argued that the total dis- 
tortion cannot be split into an elastic and plastic part. In particular in section 5 we argued that for 
the straight disclination (5.3) cannot be split because it contains the angle <p. For the dislocation 
model, however, we ignore this point of view and say that the total distortion can be split into an 
elastic and plastic part. This approach leads to Mura's new concepts of "plastic distortion" and 
"elastic distortion" of a disclination line. In particular, for the straight disclination, we assert that 
(5.3) can be split. The logical terms for the plastic part are those that are concentrated on the sur- 
face S, which correspond exactly to (8.2), or the dislocation loop density of section 4. So we find 
the elastic distortion from (3.13), (5.3), and (8.2), 
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This is the elastic distortion of the "dislocation model" of the straight disclination line. It also 
represents Mura's "elastic distortion" of the straight disclination line. Note that these expressions 
still contain the angle 0. For consistency it is straightforward to check that (8.7) and (8.4) satisfy the 
field equation (3.15). From the elastic distortion we can then find all the other elastic fields of the 
dislocation model. 

The elastic strain is obtained by substituting (8.7) into (3.16) and the result is identical to (5.5). 
Thus the elastic strain and hence stress of the dislocation model is identical to that of the corre- 
sponding discrete disclination, as we already discussed in section 119. 

The elastic rotation is found from (3.22) and (8.7) 
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In general the elastic bend-twist is given by, (113.14), 
For our special case (8.8) we find the components 
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where we have used (B3-4) and (B8). The elastic bend-twist of the dislocation model consists of the 
sum of the elastic bend-twist of the corresponding disclination line, (5.9). and the disclination loop 
density or Mura's "plastic rotation," (4.4). For consistency it can easily be verified that (8.10) and 
(8.5) satisfy the field equation (5.11), by (B8) and (Bll-12). 

Mura's new concepts therefore have a clear interpretation in terms of the dislocation model. 
Mura's "plastic and elastic distortion" of the disclination line are the plastic and elastic distortion 
of the corresponding dislocation model. Mura's "plastic rotation" is the excess elastic bend-twist 
of the dislocation model over the corresponding disclination line. 

9. The Finite Disclination Dipole 

The results for a finite dipole can be derived simply by combining the fields of a positive and 
negative basic disclination which are a finite distance apart. We choose the positive disclination 
to pass through the point (/., 0, 0) and the negative one through (— L, 0, 0), both parallel to the x?> 
axis. 

9.1. The Biaxial Dipole 
a. Geometry 

We find the defect densities for the finite straight biaxial disclination dipole from the basic 
defect densities of the discrete straight disclination line, given by (4.11-12), using the above pro- 
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cedure. The only nonvanishing components are 

cxg = n*xs[8(xi -L)- 8( Xl +L)]8 (*,) , (9.1a) 

a# 2 = - fti* 3 [8(*, -L)-8(xi + L)]8(x 2 ) 9 (9.1b) 

0g= a q [8(xi -L)-8(xi + L)] B(xt). (9.2) 

It is interesting also to find the dislocation model for this defect since it has exactly the same 
total displacement, total distortion, elastic strain, and stress. The defect densities for the disloca- 
tion model of the finite biaxial disclination dipole can be found from the basic defect densities of 
the dislocation model of the discrete straight disclination line, given by (8.4-5). Using the same pro- 
cedure again, we find 

afi = - Of[ff(*i -L)- H(xi + L)]8(x 2 ) , (9.3a) 

a§i = 0, (9.3b) 

a§i = n 2 *3[S(*i ~L) - 8(*i + L)] 8(x 2 ) , (9.3c) 

a? 2 = a^Hixi -L)- H(x 2 + £)] 8(*»), (9.3d) 

<*? 2 = 0, (9.3e) 

a§ 2 = -lli*,[8(*i-L)-8(*i + L)] 8(*«) + n 3 [#(*i - L) - #(*i + L)] 8(*i), (9.3f) 

a? 3 = 0, (9.3g) 

a§ 3 = 0, (9.3h) 

a§3 = - 0«[ff(x, - L) - ff(*, + L)] 8(x 2 ), (9.3i) 

fl& = 0. (9.4) 

The dislocation model consists of two parts: the same dislocation density as for the corresponding 
defect, given by (9.1), and a constant distribution over the surface strip between (— L, 0, 0) and 
(L, 0, 0). We have illustrated this case in figures 12-14. The dislocation model of the biaxial wedge 
disclination dipole, figure 14, was examined by Li [12], who called it a "dislocation wall." 

b. Statics 

The fields of the finite straight biaxial disclination dipole can be found from the fields of the 
corresponding basic defect, i.e., the straight disclination line, given in section 5, by the same pro- 
cedure as used above. We shall illustrate this by two examples. We first find the total displacement 
of the biaxial disclination dipole from (5.1) 

Qm fl { Xi -L)* + xl x\ x\ 1 

"• 4tt(1-v)L2 U ; (*,+L) 2 + *f (*i-L) 2 + *1 {x l + LY + x\\ 

+ ^[* ixi - L > X2) - ,p{xi+L > X2)+ m : ^){(x 1 -L)z+xr(x 1 +L)>+x*]_ 
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FIGURE 12. Dislocation model of the finite twist disclination 
dipole {biaxial and uniaxial), corresponding to figures 6 and 9. 

The dislocation densities are tt£«Oi[ff(xi — L)— H(xt+L)]&(xt) and a£= 

-a i x,[8(x l -L)-d(x i + L)]8(x 2 ). 




FIGURE 13. Dislocation model of the finite twist disclination 
dipole, corresponding to figures 7 and 10. 

The dislocation densities for the biaxial dipole are a [\ =a % 3 =~ CU[H(xi — L) 
-#(*!+£) ]«(«,) and a 3 D =n 2 X3[8(x,-/.)-8U, + /.)]8(x2). For the uniaxial 
dipole we have to add the dislocation density oc^ — — Cl 2 L[8(x\ — /.) -f8(xi + L)] 
8(jc 2 ), which is indicated by the two Burgers vectors at /> and — L. 




FIGURE 14. Disclination model of the finite wedge disclina- 
tion dipole, corresponding to figures 8 and 11. 

The dislocation density for the biaxial dipole is a D = il 3 [H(x l — L) — //(xi + 
L)]8(x 2 ). For the uniaxial dipole we have to add the dislocation density <** 2 = 
n 3 A[8(*i — L) + 8(x l + L)]8(x 2 ), which is indicated by the two large Burgers 
vectors at L and -/>. 
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where (p(x u x 2 ) is denned by (B2). Second we find the stress of the biaxial wedge disclination dipole 
from (5.12), setting (l l = [l 2 = 0, 



Gil 3 


ri (*,-L)*+*i 

.2 (*,+L) 2 +*l 


-y- 2 
*2 


*2 


11 2ir(l-v) 


( Xl -L) 2 + x 2 2 


(*i+L) 2 +*l. 


GCl 3 

<T22 ~2tt(1-v) 


ri {x,-LY+xi 

2 (xi+W+xl 


(xi-LY 
{x x -L) 2 + x\ 


(*. + £) 2 " 
( Xl + LY + xl 


G(V {x.-LY + xl 
°" 33 2tt{\-v) {x^LY + xV 






„ Gfi 3 A;2 




Xr + L I 

x,+Z,) 2 + x 2 J' 




° 12 2tt(1-* 




5^3=0, 






oft=0. 









■]• 



(9.6a) 
(9.6b) 
(9.6c) 

(9.6d) 

(9.6e) 
(9.6f) 



These results agree with the stress at large distances from a finite dislocation wall found by Li [12], 
except for a 90° rotation of the X\X 2 plane. 

Finally, to compare the finite with the infinitesimal dipole, we let the two basic disclinations 
approach each other. When L^> in (9.5) the displacement of the finite dipole approaches 
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These results agree with the displacement of the infinitesimal dipole (7.6) if we put f r = (2L, 0, 0). 
This was to be expected since this dipole is built up by confining the translation of the basic 
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disclination to the x\ axis. Next we let L 
dipole approaches 



in (9.6) and find that the stress of the finite wedge 
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These results should agree with the stress of the infinitesimal biaxial wedge disclination dipole, 
which we did not calculate in section 7. However, we noted there that the fields of such a dipole 
are the same as those of an edge dislocation. If we substitute the translation found above, £ r — (2L, 
0, 0), into (7.17) we find 6/=(0, — 2in 3 , 0). Now (9.8) agrees with (3.19) with this particular value 
for the Burgers vector. 

9.2 The Uniaxial Dipole 

We find the finite straight uniaxial disclination dipole by moving the axis of each basic 
disclination in section 9.1 back to the x 3 axis. We shall sketch out how to get results for this case. 
In general the axis of a discrete disclination line is translated from the point x" to the point x"' by 
adding to it a discrete dislocation line whose dislocation density is given by, (115.24), 



a p i(T) = 8 p (L)e, q Jl q (x r -x°'). 



(9.9) 



In our special case x°' = 0, because we move the axis to the origin. For the positive disclination 
passing through (L, 0, 0) and the negative through (— L, 0, 0) parallel to the x 3 axis we find that 
the only nonvanishing components of d p (L) are 



8 3 (L)=8(* 1 -L) 8(s 2 ), at*;>= (L, 0,0), 
d,(L) =8(*i 4- L) 8(* 2 ), atx ( ;.= (-1,0,0). 



(9.10a) 



(9.10b) 



by a coordinate translation of (A15). Hence, we find that the only nonvanishing components of the 
additional dislocations (9.9) are 



^ 2 = n 3 L[8(jci-L) + 8U 1 4-X)]8U 2 ), 
< = -n 2 L[5U,-L) +8(*! + L)]8(*2). 
We now find the dislocation density from 
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(9.11a) 
(9.11b) 

(9.12) 



Note the change in sign between this relation and (7.12). Relation (9.12) gives the dislocation density 
of the finite uniaxial disclination dipole if a^ is given by (9.1), but also of the corresponding dis- 
location model if afi is given by (9.3). The disclination density of the uniaxial dipole remains the 
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same as for the biaxial dipole, (9.2) or (9.4). The additional dislocations (9.11) are also illustrated in 
figures 13-14. 

The static fields of the uniaxial dipole can be found simply by adding the fields of the dis- 
locations (9.11) to the fields of the biaxial dipole in section 9.1. The procedure is straightforward 
and we leave the details to the reader. 

10. The Compensated Disclination Line 

As discussed in section 1110 the compensated disclination line is obtained by eliminating the 
dislocation loop density, (II10.1), 

0*=O, (10.1) 

from the defect surface S that spans the disclination line L. We shall illustrate this for our special 
case of a straight disclination. There are of course many possible choices for the surface S. We 
shall choose two of them, first the one already used and illustrated in figure 1, and second the 
semi-infinite x 2 x 3 plane for positive x 2 . 

From (10.1) and (4.5) we see that the plastic strain also vanishes, (II10.3), 

e p kl = 0. (10.2) 

With the choice for S as in figure 1 the disclination loop density remains the same as before, (4.4), 
and becomes the plastic bend-twist, (II10.4). So it has the only nonvanishing components 

<=<=n*ff(-*i)8(s 2 ). (10.3) 

The defect densities of the compensated disclination are in general given by, (II10.5), 

a pl (r) = 8 pl 8 k (S)£l k -8 l (S)£l p , (10.4) 

d pq (r) = 8 p (L)(l q . (10.5) 

For our special case we find the only nonvanishing components 

a n = -(hH (-*i)8(* 2 ), (10.6a) 

a l2 = £liH(-xi)8(x 2 ), (10.6b) 

a 32 = £l 3 H{- Xl )8(x 2 ), (10.6c) 

a 33 = -n 2 H(-x 1 )8(x 2 ), (10.6d) 

3q =£l q 8(x l )8(x 2 ), (10.7) 

where we have used (A15-16). The compensated disclination consists of the same disclination line 
along L as the corresponding discrete disclination, given by (4.12), and a constant dislocation dis- 
tribution over the surface S, which is opposite to that of the dislocation model, given by (8.4). The 
latter distribution is a simple rotation wall, as we saw. The components «i2 and a 32 represent tilt 
walls, while the components ot\\ and ot 33 represent a twist wall. 

For consistency it can be verified that (10.1-3) and (10.6-7) satisfy the definitions (4.13-14), as 
well as the relations (4.15-16). It is also easily verified that the continuity conditions (4.17-18) are 
satisfied by (10.6-7), in agreement with (II10.7-8). 

640 



The incompatibility tensor is found by substituting (10.6-7) into (4.30) and we find 

T?pg=0. (10.8) 

In general we know that the total displacement and elastic strain vanish for the compensated 
disclination line, (1110.11-12), 

4=0, (10.9) 

e k i = 0. (10.10) 

The last result (10.10) is consistent with (10.8), because the incompatibility is the source of elastic 
strain. Since the elastic strain vanishes we see by Hookes' law that the stress also vanishes. The 
only elastic field that does not vanish is the elastic bend-twist. It is in general given by, (1110.13), 

**« = -*£,, (io.il) 

and for our special case has the only non-vanishing components 

K 2q = -Cl q H(-x 1 ) 8(* 2 ), (10.12) 

by (10.3). 

We note that if we add the results of the compensated disclination to those of the dislocation 
model in section 8 we get the results for the disclination line in sections 4-5, as we already pointed 
out for the general case in section 1110. For example, we have, symbolically, (10.1) + (8.2) = (4.3), 
(10.3)+ (8.1)= (4.4), (10.6) + (8.4) = (4.11), (10.7)+ (8.5)= (4.12), (10.8) + (8.6) = (4.31), 
and (10.12)+ (8.10)= (5.9). 

We now investigate the other choice for the surface S, namely the semi-infinite plane lying in 
the x>iX\\ plane for positive %t and with a normal in the negative %\ direction. This geometry can 
easily be obtained by rotating the plane S of figure 1 clockwise by 90°. Then the only nonvanishing 
component of 8,, (S) becomes 

8i(S)=-8(*0 H(x 2 ), (10.13) 

instead of (A16). Hence from (4.2) the only nonvanishing components of the plastic bend-twist and 
disclination loop densities become 

K? q =i P r q =il q 8(x 1 )H(x 2 ), (10.14) 

instead of (10.3). Furthermore from (10.4) the only nonvanishing components of the dislocation 
density now become 

a 2 i = n 2 8(xi) H(x 2 ), (10.15a) 

£*:« = n 3 8(z,) H(xt) 9 (10.15b) 

e* 22 = -n, 8(*0 H(xi) % (10.15c) 

a3s = -fti8(*i) H(x 2 ), (10.15d) 

instead of (10.6). Finally the elastic bend-twist is also changed according to (10.11) and (10.14). 
All other relations remain the same for either choice of the surface S. 
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1 1 . The Bent Dislocation Wall 

In section 10 we found the fields of the compensated disclination for two different choices of 
the defect surface S. If we take the difference between the results for these two geometries all the 
fields vanish, except the plastic and elastic bend-twist, the disclination loop density, and the dis- 
location density. The latter becomes from (10.6) and (10.15) 



an = a 2 #(-*,)8(£ 2 ), 

a 2 i = Cl 2 S(xi) H(x- 2 ), 
a 31 = Cl 3 8(xi) H( Xi ), 

«i2 = — CliH(— x x ) 8(x 2 ), 
olti = — £l\h{x\) H(x-2.), 

«32 = — fi;)f/(— Xl) 8(x 2 ), 

ais = 0, 
a-23 = 0, 
a 33 = - flMxi) H(x 2 ) + ft 2 H(- Xl ) d(x-,). 



(11.1a) 
(11.1b) 

(11.1c) 
(Hid) 

(ll.le) 
(11. If) 

(11.1a) 

(ii. ih) 
(li.ii) 



What kind of defect is this? It represents two semi-infinite dislocation walls that meet at the Xs 
axis, i.e., it is bent dislocation wall. It has no long range stress or elastic strain. The three different 
cases are illustrated in figures 15-17, where for clarity we have shown the discrete rather than the 
constant dislocation distribution. Figures 15-16 show the geometry that has been observed by 
Schober and Baluffi [13], their figures 6 and 7. 





FIGURE 15. The bent dislocation wall with Hi rotation, 
described by the dislocation density ai 2 = — 0,iH( — Xi) 
8(x 2 ) and a22 = OL 3 3 = —Cli8(xi)H(x 2 ). 



FIGURE 16. The bent dislocation wall with Ck rotation, 
described by the dislocation density a u = 0:33= (1 2 #( ~X\) 
8(* 2 ) and a 2 i = Cl 2^ (xi)H(x 2 ). 
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FIGURE 17. The bent dislocation wall with 3 rotation, 
described by the dislocation density a 3 i = Cl^(x\) H(x-z) 
and a 32 = — n.3//( — X\)&(x 2 ). 




FIGURE 18. Edge dislocation terminating at a jog on a 
wedge disclination. 



12. Example of a Terminating Dislocation 

In a recent publication [14] we gave a concrete example of a discrete edge dislocation termi- 
nating on a jog in a discrete wedge disclination. We claimed that example illustrated the continuity 
equation (4.17). In this section we give a detailed analysis of the geometry of this example. Then we 
can verify the continuity equation for this case explicitly. 

We have already derived the general expressions for a discrete dislocation terminating on a dis- 
crete disclination. The defect densities are given by, (115.25-26), 

Mr^MLMfti + ei^ter^ 

(12.1) 



6 pq (T) = 8 p (L + L')Cl q 



(12.2) 



Here the disclination line lies along the closed curve L and /,', and the dislocation line lies along the 
connecting link L" . 

We choose the special case illustrated in figure 18. We set bi=x$'<=0, and for ft, 7 and x Q r we 
choose the only nonvanishing components ft{ andxJ] > 0. Furthermore we choose for/- a line parallel 
to the positive X3 axis passing through the point (0, x§, 0), for L ' the negative #3 axis and the small jog 
on the positive x% axis, running from the origin to £§, and for L" the rest of the positive x 2 axis beyond 
X2. Then the only nonvanishing components of the delta functions in (12.1-2) are 



83(L) = d(xi)8(x2-x°)H(x 3 ), 
82{L') = 8(x l )H(x 2 )H(x»-x 2 )8(x3), 
8 3 (L') = 8(x 1 )8(x 2 )H(-x 3 ), 
8 2 (L") = 8(x 1 )H(x2-4)8(xs). 
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(12.3a) 
(12.3b) 
(12.3c) 
(12.3d) 



Therefore we find for (12.1-2) the only nonvanishing components 

a 2 i = -0.3*2 SUO H(x 2 ) H(x? 2 -x 2 ) 8(x 3 ) -ft 3 * 2 °8Ui) H(x 2 -x? 2 ) 8(* 3 ), (12.4) 

0«s = O,8(*i) H(x 2 ) H(x*-x*) 8(* 8 ), (12.5a) 

9sa = Os8(*i) 6{x*-x*) H(x3) + n 3 8( Xl ) 8(x 2 ) H(-x 3 ). (12.5b) 

where we have used (A8). These relations represent a wedge disclination with a jog near the 
origin from which an edge dislocation emerges. Specifically, the disclination consists of a straight 
wedge disclination approaching the origin along the negative x 3 axis, (the second term in (12.5b)), 
a small jog of twist disclination running along the x 2 axis from the origin to x%, (12.5a), and a straight 
wedge disclination running parallel to the positive x 3 axis from the point x? 2 on the X2 axis, (the first 
term of (12.5b)). The dislocation consists of the following: The first term of (12.4) represents a 
dislocation that grows linearly in strength along the x 2 axis from the origin to the point x%, and 
can be regarded as an integral part of the element of discrete twist disclination. The second term 
of (12.4) is an edge dislocation running along the positive x 2 axis starting from the point x®. From 
(12.4) this edge dislocation has the Burgers vector 

bi = -£lsxl (12.6) 

Note that the effect of the terminating dislocation is to move the disclination as well as its axis 
from x 2 = for L' to x 2 = x? 2 for L. 

Now we shall show that (12.4~5) satisfy the continuity equations (4.17~18) 

OLpl s p + tlpqOpq = 0, (12.7) 

m ,p = O. (12.8) 

We need only consider the cases / = 1 in (12.7) and q — 3 in (12.8), for all other cases are identically 
satisfied. We have for (12.7) 

a 2 i,2 + 023 = -ft38(*i) H(x 2 ) H(x? 2 -x 2 ) 8(* 8 ) + ^3*2 8(*i) H{x 2 ) d(x? 2 -x 2 ) 8(x 3 ) 

-n 3 x^8( Xl ) 8(x 2 -x? 2 ) 8(* 8 ) +fl 3 8(^i) H(x 2 ) H{x\- x 2 ) 8(x 8 ) 

= 0, 

by (A6-7) and the assumption that x? 2 > 0. We have for (12.8) 

023,2 4- 033,3 = Oa8(*i) B(x t ) H(? 2 -x 2 ) B(x z ) -a 3 8( Xl ) H(x 2 ) 8(x? 2 -x 2 ) 8(x 3 ) 

+ ft 3 8(*i) 8(x 2 -x\) 8(xs) -a 3 8(xi) 8(x 2 ) 8(x 3 ) 
= 0, 

similarly. So we have shown here by a detailed analysis that a dislocation can terminate on a dis- 
clination, and that this is consistent with the continuity equations. 
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13. Summary 

We have specialized results derived for the general theory to the special case of straight 
disclination lines. First we recalled the work of Volterra and reviewed the straight dislocation to 
introduce the subject. 

We derived the geometry of straight disclination lines. We found the defect loop densities 
which are the same as Mura's "plastic distortion" and "plastic rotation." They are concentrated on 
the defect surface, which is a semi-infinite plane terminating at the disclination line. Then we 
identified the basic plastic fields (plastic strain and bend-twist) which are also concentrated on 
the defect surface. We found the defect densities (dislocation and disclination) which are singular 
along the disclination line. We noted that there is a finite dislocation density for the discrete dis- 
clination line, and we gave it an intuitively reasonable interpretation. Our results agree completely 
with Anthony's. We checked the consistency of our specific results by showing that they satisfy 
the necessary general equations: the definitions of the defect densities in terms of the basic plastic 
fields, the relations between defect loops and densities, and the continuity conditions. We illustrated 
the significance of the total Burgers vector as the displacement jump at the origin by showing that 
it vanishes for our special choice, whereas the displacement jump condition gives a result that is 
intuitively clear from the illustrations. The incompatibility tensor was found and shown to satisfy its 
continuity condition. 

Then we found the isotropic static results for straight disclinations. The displacements of a 
straight disclination were written down. Our point of departure was to regard these expressions as 
an educated guess. By differentiating we found the total distortion. From this total distortion the 
basic total fields (strain and bend-twist) can be derived. They consist of two types of terms: the first 
terms are singular on the defect surface and correspond exactly to the basic plastic fields derived 
independently above and the second terms are continuous away from the singular disclination line. 
This then provides a natural separation of the basic fields into elastic and plastic parts. In this man- 
ner we found the elastic strain and bend-twist of a straight disclination line. We showed that they 
satisfied the field equations with the defect densities. The stress followed from Hooke's law. It 
satisfied the equilibrium equations, including at the singular defect line, thus confirming that the 
original displacements were correct. 

Then we showed how the fields were modified when the disclination axis was moved away from 
the origin. The dislocation density and total Burgers vector were changed whereas the disclination 
density and Frank vector remained the same. The results can be interpreted by saying that the 
axis is moved by adding a dislocation to the disclination line. This observation facilitates finding the 
static fields of the disclination with the new axis. The fields are found by adding those for a dis- 
location to those for a disclination with its axis thru the origin. 

We derived results for the infinitesimal disclination dipole, which can be biaxial or uniaxial. 
We found the defect densities and displacements of both types and the relations between them. 
We noted that the displacement of the straight uniaxial dipole does not contain a wedge com- 
ponent. Therefore the static fields of the straight biaxial wedge disclination dipole are the same as 
those of a straight edge dislocation, confirming a result published by Eshelby. 

The dislocation model of a straight disclination was found to be a terminating dislocation wall 
(tilt and twist). We found that the dislocation model and its corresponding disclination line have the 
same incompatibility tensor, displacement, total distortion, elastic strain, and stress. The total 
distortion contains terms that are singular on the defect surface, which correspond exactly to the 
dislocation loop density derived independently before. If these terms are identified as the plastic 
distortion, as Mura did, the dislocation model follows. In this way we found the interpretation of 
Mura's "elastic distortion" of the discrete disclination as the elastic distortion of the corresponding 
dislocation model. Finally we interpreted Mura's "plastic rotation" as the difference in bend-twist 
between the disclination and its dislocation model. 

Next we found the defect densities of the finite disclination dipole. It was also possible to con- 
struct the dislocation model for this defect. The displacement of the biaxial dipole was found. It 

645 



approaches the displacement of the infinitesimal dipole when the width of the dipole approaches 
zero. The stress of the biaxial wedge disclination dipole was found to correspond to an earlier 
calculation of Li for the finite dislocation wall. 

We found the defect densities for the compensated disclination line, using two different defect 
surfaces. From these results we constructed the bent dislocation wall, which has been observed by 
Schober and Baluffi. These defects have no displacement, elastic strain, or stress. 

Finally we showed in analytical detail how a discrete edge dislocation can end on a jog in a 
discrete wedge disclination. The jog is a small segment of twist disclination line, which contains a 
finite amount of dislocation density. We showed that the analytic results satisfy the continuity 
conditions. 

In appendix A we show how the general Dirac delta functions for a curve and a surface, which 
are defined for and occur in the general theory, are specialized to the case of a straight line and a 
semi-infinite plane. It is essential to know this type of specialization in order to apply the general 
theory to particular geometries. We also calculate some useful Burgers circuit integrals. 

In appendix B we show how to take the generalized derivatives of the cylindrical coordinates. 
We carefully define the cylindrical angle <p in terms of the Cartesian Coordinates X\ andx 2 , restrict- 
ing it to the range (—77, it). Its derivative is regarded as a generalized derivative and contains a 
delta function. Also functions whose degree of homogeneity is — 1 or less will have generalized 
functions in their generalized derivatives. The usefulness of these generalized function tech- 
niques is that they allow us to carry through the traditional field operations, such as differentiation, 
even at singular lines and surfaces, in the field theoretic framework that is traditionally set up for 
continuous well-behaved functions only. This simplifies the work for discrete defects considerably, 
since we can apply straightforward mathematical techniques at the singularities, instead of intuitive 
notions as is done classically. 

In appendix C we show how the isotropic general results of III can be used to find the functional 
dependence of the static fields for straight discrete defects. 

14. Appendix A. Evaluation of the Delta Function on a Straight Line and a 

Semi-Infinite Plane 

The Dirac delta function 8(x) is defined by, (IIB1), 

r b (0, if x < a, 

8(x-x')<p(x')dx'=\(p(x), ita<x<b, (Al) 

Ja [0, ifb<x. 

This relation is illustrated by the following special cases 

§(x-x')dx' = l, (A2) 

[™ §{x-x')dx' = H{x), (A3) 

f° 8(x-x')dx' = H(-x), (A4) 

where we have introduced the Heaviside function defined by 

r// x _ f 0, if x< 0, , . . 

The derivative of the Heaviside function is easily found from (A3) by a partial integration 
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r 



„,, , BH(x) 

H (x) = — - 



dx 



= f X j-8(x-x')dx' 

Jo dx 

= — I — —. 8 (x — x ) dx 
Jo dx 



= -[8(*-*')]ff 

= 8(x), (A6) 

since 8(x) = whenever x ^ 0. From (Al) we see that we also have 

8(x-x')<p(x')=8(x-x')<p(x). (A7) 

A particular case of this relation is 

(x-x')8(x-x')=0. (A8) 

By differentiating this relation we find 

(x-x')8'(x-x')=-8(x-x f ). (A9) 

The three-dimensional Dirac delta function 8(r) is defined by, (IIB2), 

S(r — r') = 8(x\— x[)8(x> — x' 2 )8(x^ — x' 3 ). (A10) 

The Dirac delta functions for a curve L and a surface S are defined by, (IIB4-5), 

8*(L) = f 8(r-r')^L;, (All) 

8i(S) = J 8(r-r')dS{. (A12) 

We want to evaluate these quantities for the special case when L is the straight line lying along 
the X3 axis and S is the semi-infinite plane lying in the X1X3 plane for negative X\ and with a normal 
in the negative x 2 direction, as illustrated in figure 1. Then we have in particular for L 



<£L; = (0,0,dZ^) : =(0,0,db4)-, 
r'= (0,0,4), 

and for S 

dS,' = (0, dS 2 \ 0) = (0, - dx[dx! t , 0) , 
r' = (a;,', 0, x,'). 
We find that the only nonvanishing components of (All- 12) are 



(A13) 



(A14) 



8 3 (L)= [" 8(xi)8(xz)8(x 3 -x' 3 )dx^ 

J -co 
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= 8(*i)8(* a ), (A15) 

8 2 (S)= — J 8(jci — x' l )d(x 2 )d(x3 — x' 3 )dx[dx' 3 

= -H(-x 1 )8(x 2 ), , Alfi v 

where we have used (A2) and (A4). We note that these results satisfy Stokes' theorem, (IIB26), 

e ijfc 8 u (S)=-8*(L), (A17) 

by (A6). This was to be expected since L is the boundary of S. Next we want to consider some 
integrals along a curve A. that links L. We take the special case where A. is a circle of radius a in 
the X\X 2 plane centered at the origin, also illustrated in figure 1. Then we have 

x} + xg = a 2 (A18) 

on k. We wish to evaluate three typical closed line integrals on A involving generalized functions. 
The first one is 

lf(x i )8(x2)dL 1 = -(i)f{x l )8(x2)x2{a 2 -xiy^ 2 dL 2 

= 0, (A19) 

where we have used (A18) to change the variable of integration, and (A8). The second one is 

if(x 1 )8(x 2 )dL 2 = [V([a 2 -*f] 1/2 ) 8(x 2 )dx 2 + J~ a f(- [a 2 -x!YI 2 )8(x 2 )dx 2 

= f(a)-f(-a), (A20) 

where we have used (A18) and (Al). The third one is 

jj(x 1 )S , (x 2 )dL 2 = -if f (x 1 )8(x 2 )dL 1 

= 0, (A21) 

by a partial integration and (A19). Finally, for the sake of consistency we wish to show that the 
results (A15-16) satisfy (IIB15). Let cr be the circular plane bounded by A (fig. 1). Then we have 
from (A15) 

I 8 3 (L)dS 3 = I I 8(xi)8(x 2 )dxidx 2 

= 1, (A22) 

by (Al). Next we have from (A16) 

8 2 (S)dL 2 = -<b H(-x 1 )8(x 2 )dL 2 



A 



=-H(-a) + H(a) 

= 1, (A23) 



by (A20) and (A5). So these results agree with (IIB15). 
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15. Appendix B. Generalized Derivatives of Cylindrical Coordinates 

The expressions that are differentiated in the text contain the cylindrical coordinates p and 
<p. The geometric relation of these coordinates to the Cartesian Coordinates X\ and x± is illus- 
trated in figure 5. The functional relation for p is 

p(*i,**) = (*J + sj) 1 fr. (Bl) 

The relation for (p is customarily given as tan -1 (#2/^1 )• However, this leads to some ambiguity for 
two reasons: (1) The arctangent is usually taken to be restricted to its principal value, lying in the 
range (— 7r/2, 7t/2), especially in integral formulas. (2) The angle ip is usually regarded to be multiple- 
valued, a procedure that conflicts with the assumption of uniqueness for the functions that occur 
in the general defect theory. Therefore we shall arbitrarily restrict <p to the range (—77, 77), so 
that it is uniquely defined everywhere, but possesses a discontinuity of 27T across the negative X\ 
axis, as shown in figure 5. This choice is convenient for the calculations in the text, for then the 
discontinuity occurs at the semi-infinite plane S of appendix A (fig. 1). So the functional relation 
for (p is 

ip(x u x 2 ) = tan- 1 {x 2 lxi) + 7rH(—Xi)[H(x2)—H(—x 2 )] 9 (B2) 

where the arctangent is restricted to its principal value and the Heaviside function H is defined 
by (A5). We see that (p, so defined, assumes the values indicated in figure 5. 

In order to maintain the necessary uniqueness when differentiating the discontinuous function 
<p, it becomes necessary to take the generalized derivative. The discontinuity will then lead to a 
generalized function, the Dirac delta function. We refer to Gel'fand and Shilov [8], chapter 1, 
section 2.2, for the detailed procedure on taking this generalized derivative. The result is 

<P,i = -Wp 2 , (B3) 

<Pfi = xilp 2 +2irH(-xi)8(x 2 ). (B4) 

Intuitively, these results could have been anticipated. They contain the ordinary derivatives of 
the arctangent function plus a term having a delta function across the negative X\ axis, which 
originates from the discontinuity of <p. 

By contrast the generalized derivatives of p and In p are the same as the straightforward 
classical ones. For the latter we have 

(lnp),i = Wp 2 , (B5) 

(lnp), 2 = Wp 2 . (B6) 

Another function whose generalized and classical derivatives coincide is the following 

\ P ) ,y P P 

where we have used the convention that Greek subscripts have the range 1, 2. When written out 
in full these equations are 

Wj P 4 ' V/,2 p 4 ' 

\ p* /,, p 4 p 2 v p 2 ; >2 p 4 p 2 
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VpV,. p 4 ' \p 2 /,2 V "' 



The above derivatives (B5— 7) are all locally integrable in the X\X 2 plane, including the singularity 
at the origin. This is the reason why the generalized and classical functions coincide. 

However, if we take further derivatives of (B7), the results (having a p~ 2 dependence) will no 
longer be locally integrable at the origin, because then the algebraic singularity has the same 
strength as the dimension of the space (i.e., 2). To maintain uniqueness it becomes again necessary 
to take the generalized derivative. We refer also to Gel'fand and Shilov [8], chapter III, section 3.3, 
for the detailed procedure on taking this type of generalized derivative. The results will contain 
expressions homogeneous of degree — 2 (such as p~ 2 ), which have to be regarded as generalized 
functions. One result is 

(^)=^-2^+-8^(p), (B8) 

where we have defined the two-dimensional Dirac delta function 

8(p)=8(*i)8(*«). (B9) 

When (B8) is written out in full we have 

Xi \ 1 x\ l ocx \ #1*2 



X2 

7 2 / 



p 4 \P 2 /,2 P 2 P 4 



Another result is 



V p 4 7,6 p 6 p 4 4 



or, when written out in full 



=-4^- 2 
P 4 /,2 p 6 ' 



/ X 3 t\ x\ n x\ 3tT , ^ 

\ P 4 7,i P 6 P 4 V P 4 7,2 p 6 P 4 4 vp " 

\ p 4 7,i p 6 p 4 4 KH ' \ p 4 7 >2 p 6 p 4 
\p 4 7,i p 6 \p 4 7,2 p 6 p 4 4 

Further generalized derivatives can be taken similarly, leading to homogeneous generalized 
functions of degree — 3. Two results we need are 

{?), a = - 2 p-« 8Aph (B11) 
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\ p 4 A 



= -4 



XgXpXy xJJsy + Xp&yg 77 



P 4 A* P° 

The latter equation written out in full is 

3tt 



(8 n f$Mp)+Sev&.a(p) + 8y<£.fi(p)). (B12) 






X\X\ X\ 7T 



V P 4 Ai P 6 p 4 4 ' VK/ 

/ ac« \ X\X\ TT 

(j0„=- 4 ^-4 8 "«»- 



— *-^ + ^-4 8 ' ( ' ,K 



=-4-f+2— — S, 2 (p). 
,2 p b p 4 4 



16. Appendix C. Derivation of Static Fields of a Straight Defect from the 

General Theory 

In III we derived general expressions for the static fields in an infinitely extended isotropic 
body. We showed there that a condition on the plastic strain for the integrals to converge is that 
e£(r) approaches zero faster than r _1 as r^°°. This condition is not satisfied for straight defect 

lines, as can be seen from (3.3) and (4.7). However, it is still possible to find the functional depend- 
ence of the static fields from the general formulas by ignoring infinite terms that appear in the in- 
tegrations. The correctness of the resulting static fields can then be assessed by checking them in 
the general field theory, as is done in the text. 

The general expressions that we shall use are the following: For a dislocation line, the dis- 
placement, (1114.21), 



af(r)=— f iRMbkdS'i+RjkkbidS' — R Mk bkdS' 3 + - !L -Rjwb* 

SV J s L 1 — v \—v 



is: 



(Cl) 



and the elastic distortion, (1114.23), 

bi 



/Mr) = 



877 J L 



R t lnn€ijkdL k + R jcnnGikldLj + 



R 



P 



&ijR,knn I €klmdL' m 



For a disclination line with its axis through the origin, the displacement, (1116.10), 



uf(r) =— €kmn£lmX' n 

07T Js 



1 v 

R jjk dSj + R jiiSikdS '. — R , ijkdS! + R jjjdS L 



l-v 



\-v 



(C2) 



(C3) 



and the basic elastic fields, i.e., the elastic strain, (1116.11), 

etj (r) = q— <j) €ipqfilpXq ' R ,knn(^ikldLj )(,j) + K _^ dijR ,knn \^klmdL' m \ + R , nn {£ljdL\ ) (jj) 

+ \T~Z ^y^ ' wn \£l m dL' m , 

and the elastic bend-twist, (1116.14), 

Klj (r) = — <b [ € kpqQ pXq {R iklnndL j — f R ,ji nn dL JJ — R ,knn€iklfydL 'J , 



(C4) 



(C5) 
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where 

R=(X i X i yi\ (C6) 

X; = Xi-x\. (C7) 

The expressions (Cl-5) are line and surface integrals. For the special geometry of figure 1 the 
curve L is given by (A13) 

dL\= (0, 0, dI4) = (0, 0, d4), r'= (0, 0, 4). (C8) 

Hence we have 

R={p*+Xiy'\ (C9) 

P ={x\+x>yi\ (cio) 

by (Bl). The surface S is given by (A14) 

dS[ = (0, dS! 2 , 0) = (0, - dx[ dx' z , 0) , r' = (x[ , 0, *£ ) . (Cll) 

Hence 

R=(JfJ + x|+^)^. (C12) 

We shall first investigate the line integrals. All the necessary results can in general be derived 
from differentiating the basic integrals / RdL[ and / x\ RdLj . From (C8) the only non-vanishing 
components of these integrals are / RdL s and / x^ RdL% . Since these integrals are along an in- 
finite straight line they will diverge. So we first find the finite integral 



Rdx' 3 = RdXs 

J-L+X 3 J-L 

= [ipZ\n(X 3 + R)+lx 3 R]'l L 



= p* In ^ Kp ^ L ' +L (p« + L*) ^ . (C13) 

The integral we want is obtained when L — > oo, but then the terms in (C13) containing L will also 
approach infinity. However, these terms can be regarded as constants in the subsequent differentia- 
tions. Therefore they do not contribute to the fields and we can ignore them for practical purposes. 
So the functional dependence of the desired integral goes as 



I 



RdL^->q, (C14) 

L 

where we have defined the function 

q = -p 2 \np. (C15) 

The justification for this procedure can also be given as follows: A sufficiently high derivative of 
(C13) will converge for L— > °o to a finite function of p. The same derivative of (C15) will give the 
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same finite function. Now the relations (Cl-5) contain only derivatives of integrals such as/ RdL}. 
Therefore we can use (C14— 15) in them directly. We find the functional dependence of the other 
integral in the same way. First we find the finite integral 



x :i R dx g = I ( x ,i — X 3 ) R dX fi 
J-L+x 3 *-L 



= \ix,p* In {Xi + R) +ix3XJl-iR* 



and hence, comparing with (C13). 






'Rdx£ -^ x-jq. 



(C16) 



(C17) 



We now consider the surface integrals. All the necessary results can in general be obtained from 
f RdS'i and j x[RdS\. From (Cll) the only nonvanishing components are fRdS 2 , J x\RdS' 2 , and 
/ x' z RdS' 2 . The derivation of the first integral may be indicated as follows 

\ RdS: 2 ->- f I Rdx[dx' 2 ^> I {X* + x*) In (X^x*) 1 *^-* \ (X\ + x\)\n(X\+x\) 1 f % dXi 

J s J-L J-L J-L Ja?j 



{hX\ + X iX D In (Ar« + *»)>#-|^ tan" 1 - 1 -(!*? + §*,*§) (C18) 



X\ 9 



Here the first result follows from the geometry of the surface S in figure 1, the second by integra- 
tion over x' z which goes through the same as above but with R given by (C12) instead of (C9), and 
the third by a change of variable. So the functional dependence becomes 



f RdS' t -*g, 

where we have defined the function 

g = — (i*f + *i*|) lnp-h f*l^+Q*?+t*i*f). 

The next integral goes as 

f x[RdS' 2 -+ [° x[(X\ + xl) In (X\ + x\y*dx[ 

-> [ {xi-Xi)(X\+s%) In (X\ + x\)V*dXi. 

J x x 

This leads to the functional dependence 



i 



x'.RdS'^h, 



(C19) 



(C20) 



(C21) 



(C22) 
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where we have defined the function 

h = Xl g+\p 4 lnp+&x*-ix*xl. (C23) 

The last integral goes as 

\ x^RdS't^>-\ i x' s Rdx' 1 dx' 3 -+x 3 g. (C24) 

JS J-Lj-L 

The functions q, g, and h depend only on %\ and x-z. The derivatives of </that we shall need are the 
following 

<?,n=-2 lnp-1-2^, (C25a) 

P 2 

~ x,x 2 
9,12= -2 — , (C25b) 



P 2 






<7, 22 =-2 lnp-1-2 — (C25c) 



g, Qa =-4 (lnp+1), (C25d) 

(C25e) 

(C25f) 
(C25g) 
(C25h) 
(C25i) 
(C25j) 
(C25k) 



9,111 = 


-2 


X\ X\X 2 

7 2 ~ p 1 " 


9,112 = 


-2 


*2 + 4 fl*2 

P 2 P 4 


9,122 = 


-2 


r 4 

P 2 P 4 


9,222 = 


-2 


%2 . X\X 2 

4 — — 

P 2 P 4 


Q,\aa 


-4 


P 2 ' 


9,2aa 


-4 


x 2 

~9.' 



g,n««=-4 (— -2^+7rS(p)), 



„ *l£ 2 

9,i2« Q = 8— -, (C251) 



P 



/ 1 x 2 \ 

9,22aa = -4(^— -2- 2 ; +7r5(p)j. ( C 25m) 



Here the Greek indices range only over 1 and 2. The results (C25k-m) follow from (B8). The other 
derivatives we shall need are 

S,ii2 =-2 — , (C26a) 
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#,i22=-2(lnp + i+%), 






g f laa = — 4(lnp + l), 

g,2aa=4<£>, 

A ,122 = — 2*i In p +4>x 2 <p—%xu 
A ,222 = 6*2 In p + 4*i<pH-6*2-f i* 2 , 
A,iaa = — 4*i In p-f4jc 2 ^-ii, 
A,2aa = 4* 2 In p + 4*i (£> + 6* 2 — i *2 . 
The displacement of the dislocation (Cl) now reduces to 



u{ 



J_ 

877 



g,aapbf3&j2 + g,2aabi— " g,i2ab a + - g,iaab2 

Y — V \—V 



(C26b) 

(C26c) 
(C26d) 
(C26e) 

(C27a) 
(C27b) 
(C27c) 
(C27d) 

(C28) 



The individual components are 

6i 



ll Stt 



Stt 



r i i 6 2 r i 2-1/ i 

jz^^22 -*4«j - ^ [t=^ - 1^*4 



b. 



U Z ^gZcta, 



(C29a) 
(C29b) 
(C29c) 



which give (3.12) by (C26), except for some trivial constants. The elastic distortion of the dislocation 
(C2) reduces to 



/3«0 = 



8tt 









(C30a) 

(C30b) 
(C30c) 



where e a /s is the two dimensional permutation symbol, i.e., €n = €22 = 0, e V 2 = ~ 621 — 1. The non- 
vanishing components of (C30) are 



)8„=- 



877 



*■— it 



9,112 1 , 62 f g.lll 

" <?,i22 . , , 62 { q,M n 1 



(C31a) 
(C31b) 
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___b±_ \ q,vz2 
Stt 



r <?,i22 1,62 r <?,ii2 , 



o _ J&i r <7,222 1 , i 2 \ q,l22 



/3 — ^3 
P13— Q^' 2 ^' 



/3 — ^ 3 

P23— jto8> lvV9 



which give (3.14) by (C25). The displacement of the disclination (C3) reduces to 



u 1 



fia 



h ;iaa §i2 "~" h ,/22 "I" ~ 7 h ,mo 



+ 



1»" ,<-^ ■ i 
— v —v 






Stt 



x 3 g, a yYdi2 + xng , 2 yy5ia ~TZT M.fia+'r" (*3g) ,tvvS a 



The individual components are 

T ftl^3 



i v i n 2 * 3 r i ] 

T n^af i 2-* i n 2 * 3 r 1 l 

~87r[l-^' 222 "l-^' 2a 



8tt 



(C31c) 
(C31d) 
(C31e) 
(C31f) 



(C32) 



(C33a) 



(C33b) 
(C33c) 



which give (5.1) by (C26— 27), except for some trivial linear terms. The elastic strain of the disclina- 
tion (C4) reduces to 






e 3 3 = 0. 
The individual components are 



(C34a) 

(C34b) 

(C34a) 



en = — 



Qi*3 

8tt 



g411 

1-v 



"2,11' 



^2^3 r £412 

8tt 



r qm2 n 1 1 Q3 r <?,n /roc 
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e>>2 = — 



Qlg3 

8tt 



g.122 



■q,n 



8tt 



e33=0, 



_ fijXn q,ii2 to'iXa q,\ii .ftg 7, 12 
ei2 " 8tt 1-i> 8tt l-i; 8tt1-i/' 



ft-«*-] + St[w-^]' <C35b) 

(C35c) 
(C35d) 



111 ^,12 H 2 



q^ 

\-v 



q.w 



e 3 i = — 



8tt 



\-v 



q,w 



H 2 g,i2 

877 1-1/ ' 



(E35e) 



(C35f) 



which give (5.5) by (C25), except for some trivial constant terms. The elastic bend-twist of the 
disclination (C5) reduces to 



fly 

Kay— W~ taQQfivv, 



(C36a) 



Ka3 — 7T~ [ClyX3G0yq,tx0in> + £laGapq,0vv], 



K3y = 0. 



The individual components are 



K\\ — o q&w, 



K21 — ^ q,2vi"> 



K31 =0, 



_n 2 

K12 — « q:iw> 



fi 2 



*32=0, 



Kl3 = 



K23 : 



(i\xa , n 2 ^3 , n 3 

^ 9,i2t ^ H~ ^ q,iiw~r 71 9*2^9 



8tt 



q, ti w 



Q2X3 

8tt 



■9.12* 



n_ 3 



(C36b) 
(C36c) 

(C37a) 

(C37b) 
(C37c) 
(C37d) 

(C37e) 
(C37f) 

(C37g) 
(C37h) 



K33 : 



8tt 



<7,2i> 



a 2 

to* 4 ' 



which give (5.9) by (C25). 



(C37i) 
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